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ABSTRACT 

Non-Gaussian distributions of cosmic microwave background (CMB) anisotropies have 
been proposed to reconcile the discrepancies between different experiments at half-degree 
scales (Coulson et al. 1994). Each experiment probes a different part of the sky, furthermore, 
sky coverage is very small, hence the sample variance of each experiment can be large, 
especially when the sky signal is non-Gaussian. We model the degree-scale CMB sky as 
a Xn fi e ld with n-degrees of freedom and show that the sample variance is enhanced over 
that of a Gaussian distribution by a factor of (n + 6)/n. The sample variance for different 
experiments are calculated, both for Gaussian and non-Gaussian distributions. We also show 
that if the distribution is highly non-Gaussian (n < 4) at half-degree scales, then the non- 
Gaussian signature of the CMB could be detected in the FIRS map, though probably not in 
the COBE map. 

Subject headings: cosmic microwave background - cosmology: theory 



Recently, several groups have reported results on the measurement of anisotropies of the 
cosmic microwave background (CMB) at degree scales (de Bernardis et al. 1994; Cheng et 
al. 1993; Gunderson et al. 1993; Meinhold et. al 1993; Schuster et al. 1993; Tucker et al. 
1993; Wollack et al. 1993). The beam size, beam throw, most sensitive angular scale, sky 
coverage and quoted rms temperature anisotropies are summarized in Table 1. The results 
from these experiments do not agree with each other, in particular the results from the same 
experiment, MAX, for two different part of the sky, the Gamma Ursa Minor (GUM) region 
and the mu-Pegasi (MuP) region, contradict each other at 2a level. A way to reconcile these 
measurements is to have a non-Gaussian distribution of temperature anisotropies at half- 
degree scales (e.g. Coulson et al. 1994). At present, there are still large uncertainties in all 
experiments due to foreground subtractions, therefore the need to invoke non-Gaussian tem- 
perature distributions remains to be established. Since different experiments probe different 
part of the sky and the sky coverage of each experiment is small, the sample variance of each 
experiment can be large, especially when the sky signal is non-Gaussian. The goal of this 
paper is to quantify the difference in the expected sample variance between non-Gaussian 
and Gaussian fields in order to determine if this effect could be responsible for the discrep- 
ancy between experiments. We also estimate the minimum sample size for each experiment 
in order for the sample variance to be less than 20/ifT, both for Gaussian and non-Gaussian 
distributions. 

Regarding the statistical analysis, the most important quantity is the number of indepen- 
dent measurements in a single sample. However, the data points in CMB experiments are 
correlated and therefore contain less statistical information. Thus, it is useful to determine 
the effective number of data points, N e , defined as the number of independent measure- 
ments of temperature anisotropies in each experiment (we will discuss how to estimate N e 
for each experiment later). Expressed in terms of N e , the sample-averaged rms temperature 
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anisotropy is 

where ATj/T is an independent measurement of temperature anisotropy at an angular scale 
corresponding to the beam size. The ensemble average of ATj/T vanishes and the ensemble 
average of A is of ft , where a t h is the theoretical prediction for the temperature anisotropy at 
the angular scale probed by the experiment. 

Due to the stochastic nature of A will be a random quantity which can vary between 
different experiments. The sample variance is the fluctuation around af ft , the ensemble- 
average of A. In the case where is Gaussian, the sample variance is given by 

vim = (A 2 ) - (A) 2 = J-(A) 2 , (A) = al. (2) 

Here the angle-brackets (...) denote an ensemble average. The sample variance will be 
reduced by increasing the sky coverage by an amount which will scale with the effective 
number of data points N e as a sam oc \/\fW e . 

We would expect the sample variance to be larger when the temperature anisotropies 
are non-Gaussianly distributed. Since the functional space of non-Gaussian distributions is 
unlimited, one is faced with the question of choosing appropriate distributions to model the 
sky. As we have stressed before (Luo 1994), a xi distribution with n-degrees of freedom is one 
of the simplest and most natural choices. By varying n, it provides a family of distributions 
that range from highly non-Gaussian (small n) to nearly Gaussian (large n). Furthermore, 
it provides a good fit to the statistics of temperature fluctuations from global topological 
defects and non-topological defects in the framework of the 0(iV)<7-model, where a global 
symmetry O(N) is broken to 0(N — 1) by an N-component real scalar field = (0i, <f)jy) 
in the early universe (Turok & Spergel 1991). In this paper, we will model the degree scale 
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CMB sky as a x\ field, with 

AT 



E(tf-" 2 ), (3) 



T 

1 3=1 

where = l...n, are n-independent Gaussian variables with zero mean and variance a 2 . 
The ensemble average of A is related to a through (A) = 2na 4 . To calculate the sample 
variance of the rms temperature anisotropy of a x 2 field, we will utilise some of the higher 
moments of a Gaussian, i.e., 

(S A ) = 3(5 2 ) 2 , (5 6 ) = 15(5 2 ) 3 , (5 s ) = 105(<5 2 ) 4 , (4) 

and the following identities for Gaussian variables: 

n 

(J2 5 2 5 2 )=n(n + 2)a\ (5) 

ij=l 

n 

(E 5 2 5 2 5 2 k )=n{n + 2){n + A)a\ (6) 

ijk=l 



and 



( E W*i 2 > = n(n + 2)(n + 4)(n + 6)a 8 . (7) 

ijkl=l 



After some algebra we have 



/ 2 \ n + 6 2 

^sam - J —\\-jT)sam) {°) 



This is the main result of the paper. For a x\ distribution, the sample variance is enhanced 
by a factor of (n + 6) jn relative to a Gaussian distribution. As we expected, as n becomes 
much larger than 6, the enhancement is negligible. 

The effective number of data points, iV e , depends on the detailed sampling scheme. If the 
experimental data are sparsely sampled, i.e. the distance between data points is much larger 
than the beam size, then N e is approximately the number of experimental data points. If the 
data is over-sampled so that correlations between the data points are important, then N e can 
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be estimated as follows. The sample variance, a sam , of a Gaussian-distributed temperature 
anisotropy can also be expressed as (Scott et al 1994): 

CTsam = ^J dQ 1 dn 2 C 2 (9 1 , 9 2 ), (9) 

here Q is the solid angle covered by the experiment and C(9i, 82) is the two-point temperature 
correlation function. Combining equations (9) and (2) gives an expression for N e : 

N = (10) 

The data analysis of all experiments involves using the Gaussian auto-correlation function 
(GACF), where the correlation among the data points is approximated as 

C(9)=C(0)eM-0 2 /29 2 c ), (11) 

with 9 the angular distance between two data points and 9 C the coherence angle. Since the 
two point correlation function depends only on the difference between two angular directions, 
after changing variable equation fllOD can be reduced to 

N e = QC 2 {0)/ J dVLC 2 {9). (12) 

For a single one dimensional temperature scan, it is easy to evaluate N e since 

5^/^=^, iHH (13) 

where 9 m is the maximum angular difference among data points with respect to the observer. 
The solid angle covered by the experiment is Q = 9 m x 9fwhm, hence one can express N e 
in terms of Q as 

N < = T^ir-^r)- (14) 

Although the above formula is derived for a single one dimensional scan, it is also applicable 
to multiple scan experiments when the correlation among different scans is small. Estimates 
of N e for the different experiments are listed in Table 2. 



The sample variance, <J sam , is directly proportional to the theoretical prediction of the 
temperature anisotropy. The theoretical model which we used to calculate the numerical 
value of o sam is the standard CDM (Q — 1) model, with a scale invariant primordial power 
spectrum. We take the Hubble constant to be h = 0.5 and the baryonic faction fi& = 0.06, 
for consistency with the big-bang nucleosynthesis bound (Walker et al. 1991). Theoretical 
values of the temperature anisotropy at degree scales for the experiments are taken from 
White et al. (1994). 

A question of practical interest is just how large should a sample be in order for the 
experimental result to be informative? One criteria could be that the sample variance is much 
smaller the the instrumental sensitivity in AT. From equations (2) and (8), we estimate the 
minimum sample size N min that gives rise to a sample variance of 20fiK in AT, both for 
Gaussian and for x 2 distributions. Minimum sample size for current experiments are shown 
in Table 2. Sky coverages of ARGO (de Bernardis et al. 1994) and Saskatoon (Wollack et al. 
1994) is already large enough so that the sample variance is subdominant in each experiment 
if the sky distribution is Gaussian. However, substantial sky coverage is needed if the sky 
distribution is highly non-Gaussian, especially for those experiments (MAX & MSAM) that 
samples around the Doppler peak of the radiation power spectra, because the theoretical 
predictions of the temperature is much higher. 

If the CMB sky distribution is highly non-Gaussian on degree scales, the signature of 
non-Gaussianity may well be detected in the large angular scale CMB maps (Smoot et al. 
1992; Ganga et al. 1993). We will use the skewness, 



to characterize the lowest order deviation of CMB distribution from a Gaussian. The skew- 



ness for the Xn distribution is fis = y 8/n. If the sky distribution is a x\ fi^d on an angular 

scale 9 ~ 1°, after smoothing with beam size 9 S , the distribution is x\'i with n! = n(9 s /9 ) 2 . 



((AT) 3 ) 



(15) 



1^3 = 



((AT) 2 ) 3 / 2 
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For COBE, 9 S = 7° and n' « 50n, and for FIRS, 6 S = 3.8° and n' « 14n. Thus for models 
with n > 4, the skewness of temperature fluctuation on COBE and FIRS scales is at least 
li 3 > 0.1 for COBEand /i 3 > 0.30 for FIRS. The cosmic variance of the skewness /i 3 
is fx 3 ~ 0.18 for large angular scale experiments like COBE or FIRS (Srednicki 1993). For 
COBE, the cosmic variance is larger than the non-Gaussian signals, thus it will be impossible 
for the non-Gaussian signal to be detected. However, for FIRS, the smooth non-Gaussian 
field will still stand above the cosmic variance. 

Finally, we shall comment on reionization. Early reionization is likely if the density 
fluctuations are non-Gaussian and objects can form at an earlier epoch. If the universe was 
reionized at an early epoch (z ~ 100), then the degree scale temperature anisotropies are 
dramatically reduced. In this case the sample variance will also be reduced dramatically. 
As an example, in Table 2 we will also list sample variance (for different experiments) for a 
re-ionized CDM model with optical depth r = 1 (Kamionkowski et al. 1994). We conclude 
that we cannot reconcile the high detection of MAX/GUM experiment with theoretical 
predictions for any x 2 distribution. If the MAX/GUM result is confirmed, models with early 
reionization will be ruled out regardless of the statistics of CMB at half-degree scales. 

We want to thank D. Scott and B. Moore for their helpful comments. This work was 
supported by a grant from the NSF. 
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Table Caption 



Table 1: A brief summary of the experimental situation at degree scales. 
Table 2: The sample variance ^a sam of rms temperature anisotropies in various experiments, 
in units of /iK, for Gaussian and xl, sky distributions. The theoretical predictions for various 
experiments are given for a CDM dominated universe with Q = 1, Qb = 0.06, h = 0.5 and a 
COBE normalized Harrison-Zel'dovich spectrum. The numbers in bold face are for standard 
recombination and those in brackets are for a fully reionized universe with optical depth 
r = 1. 
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Table 1 



Experiment 


v r VV ll Ivi 


n 

^ chop 


Or 

v c 


sample size 


AT(fiK)(±2a ) 


ARGO 


52' 


1.8° 


30' 


63 


41 - 71 


MAX (GUM) 


0.5° 


1.3° 


25' 


165 


85 - 162 


MAX (MuP) 








21 


< 68 


MSAM (single) 


0.5° 


40' 


.5° 


14 


16 - 60 


MS AM (double) 






.3° 




30 - 85 


Saskatoon 


1.44° 


2.45° 


1.44° 


48 


28 - 57 (±l(j) 


SP91 


1.5° 


2.95° 


1.5° 


13 


13 - 45 


White dish 


12' 


23.6' 


0.15° 


5 


< 63 



Table 2 



Experiment 


(AT) th (fiK) 


N e 


mm 


N* 2 

- 1 * mm 


Gaussian 


Xn=l 


ARGO 


44(27) 


63 


47(7) 


289(49) 


19(11) 


31(18) 


MAX (GUM) 


71(38) 


14 


318(26) 


2226(182) 


44(23) 


71(38) 


MAX (MuP) 
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45(24) 


74(39) 


MSAM (single) 


76(44) 


8 


417(47) 


2919(329) 


54(31) 


88(50) 


MSAM (double) 


46(22) 


13 


56(3) 


392(21) 


29(14) 


47(23) 


Saskatoon 


38(27) 


27 


26(7) 


182(49) 


20(14) 


33(23) 


SP91 


38(27) 


6 


26(7) 


182(49) 


29(21) 


47(33) 


White dish 


41(19) 


5 


35(2) 


245(14) 


33(15) 


54(25) 
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